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Summary 


In-depth gap heating ratios, q(z)/qref» were predicted down RSI tile sidewalls 
based on temperature measurements obtained from the JSC arc-jet Wing Glove model 
tests. The object! ,'es of the study were to develop gap heating ratios which 
resulted in the best possible fit of test data and to produce a set of engineering 
verification heating ratios similar in shape to one anotner which could be used 
at various body points on the Orbiter during reentry. The Rockwell TPS Multi- 
dimensional Heat Conduction Program (XF0031) was used to perform 3-0 thermal 
analyses using a 3.0 in. thick sectio;i of a curved RSI tile with 283 nodal points. 

The results of the correlation with test data showed that the predicted heating ratios 
were significantly higher down in the gap than the zero pressure gradient values for 
T/C stacks 39 and 38 on the Wing Glove model. For stack 37 (in a low pressure 
region), the baseline heating ratio overpredicted the temperature data. This 
analysis, which showed that the heating ratios were a strong function of the 
product of pressure and pressure gradient, will be used to compare with recent 
Gap/Step and Ames Double Wedge test/analysis results in the effort to identify 
the Orbiter gap response to High AP flight environment. 

(NASA-CR-160146) PREDICTION OF IN-DEPTH GAP N79-20068 

HEATING RATIOS FPOH WING GLOVE MODEL TEST 
DATA (Rockwell International Coro., Downey, 

Calif.) 57 p HC A04/MF A01 ‘ CSCI 20D Dnclas 

G3/02 18.)31 
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1.0 Introduction 

A series of three test programs (References 1-3) was performed at the NASA/ JSC 
lOMW Arc Jet Facility using a model of the Orbiter Wing Glove region. This 
model was located some 25 In. downstream of a 40 In. diameter conical nozzle at 
the facility. During the third test program of the series (Reference 3), 
In-depth temperature measurements were obtained at a number of stacks of thermo- 
couples located In the gaps of the model. A primary result of this test program 
showed that when filler was Installed In the longitudinal (stagnation line) 
gap the temperatures In the circumferential gap increased some 1 00-40 above 
those for tests with no gap fillers. As a result, the data from the t^ve 
filled stagnation line gap tests (C-188, 189, 190, 194, and 195) was used to 
correlate with thermal math model predictions. Tests C-188, 189, and 195 
were high pressure tests. Test C-194 was a moderate pressure test and Test C-190 
was a low pressure test as seen In the surface pressure plot of Figure 1. Data 
from three thermocouple stack locations (39, 38, and 37) as seen In Reference 3 
was used for each of the five test cases. The principal surface data Including 
pressure, pressure gradient, and temperature near each of these stacks and the 
test conditions for all five cases Is presented In Table I. 

The TPS Multidimensional Heat Conduction Program (XF0031) described In 
Reference 4 was used to develop an analytical math model of one quarter of an 
RSI tile. Both the derived value of heating ratio, q(z)/qref» and the original 
zero pressure gradient ( baseline) value obtained from Ref. 5 were used In the 
analysis. Reference heating rates and surface pressures at the top of each 
stack were input as a function of time for each arc-jet test. Engineering 
verification heating ratios were developed that yielded smoothed curves, roughly 
parallel to each other, that were a function of the product of surface pressure 
and pressure gradient. The remainder of this document presents the details of 
the analytical model developed and the correlation with the Wing Glove model 
test data. 

2.0 Analytical Model 

2.1 Basic Description of Model 

A 283 node 3-D model of a one-quarter curved RSI tile as shown in 
Figure 2 was developed and analyzed with the Rockwell XF0031 thermal 
analyzer program. With this program the solution to the 3-D heat 
conduction partial differential equation Is obtained by the finite 
difference method. The method used to integrate the heat flow 
equation for all of the cases analyzed was the Crank-Nicolson method, 
which is based on a "mid" difference implicit numerical analysis 
technique. 

The model used In the analysis had an LUOO RSI depth of 3.0 In. plus 
a_0.015 In. coating, 0.16 In. SIP thickness, and aluminum thickness 
(t) of 0.063 in. The properties of these four materials (RSI, coating. 
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SIP, and A1 2024) were Input to the analytical model as either constant, 
a function of temperature only, or as a function of both temperature 
and pressure, as obtained from Ref. 5. The bottom side of the aluminum 
structure was driven at a temperature corresponding to that measured 
for each test. It was assumed the Sllfrax tile across the circumferential 
gap from the RSI tile had the same properties as the RSI tile. 


In the analytical model the specific heat ai./ thermal conductivity of the 
0042 coating and A1 2024 structure were Input as a function of temperature 
only. The RSI specific heat was also Input as a function of temperature 
only while the RSI and SIP thermal conductivity were input as a function 
of pressure and temperature. The surface pressures at each stack were 
Input as a function of test time. The SIP specific heat, coating emissivlty, 
and the density of all materials were Input as constant values. 

The radiation view factors down the gap were Input from calculations 
obtained from the 2-D Cross-String Method Program. This program was 
written by Rockwell for use on the Hewlett-Packard 9830A Mini -Computer. 

The scrIpt-F radiant Interchange factors across each enclosure were 
computed by the XF0031 Program by solving a system of linear algebraic 
equations, based on the emissivlty, nodal surface areas and Input view 
factors . 

The temperatures and materials on either side of the gaps were assumed 
equal to each other at corresponding depths (thermal symmetry). The 
surface emissivlty of all surface nodes was assumed to be 0.85. The Initial 
temperature of all nodes was taken to be 80®F. This temperature was also 
assumed to be the sink temperature Inside the vacuum chamber for all surface 
nodes radiating to the sink. The most sensitive parameters In the analysis 
were the heating ratio down the walls and the reference heating rate at each 
T/C stack investigated. These Items are discussed In detail in the next 
section. 

Z.Z Heating Ratio and Reference Heating Rate 

The Wing Glove region upon w.hich this thermal analysis was performed included 
the distance along the circumferential gap from thermocouple stacks 43, 39, 
38, 37, to 34 as seen in Reference 3. The main emphasis of the analysis 
was upon stacks 39, 38, and 37. Stacks 43 and 34 were at the junction of 
the circumferential and longitudinal gaps. Stack 43 was in a stagnation 
region, and stack 34 was in a separated flow region. 


As noted in Figure 2, two specific models of the T/C stacks were used depending 
on which stack was being Investigated: 1) Model 38-39-43 or 2) Model 38-37-34. 

The stack 39 analysis used the first model, the stack 37 analysis used the 
second model, and the stack 38 analysis used both models. When the analysis 
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was performed for stacks 38 or 39 the reference heating rate at stack 38 was 
Input at nodes 11, 12 and 13 (Refer to Figure 2), the reference heating 
rate at stack 39 was Input at nodes 6 , 7, and 8 and the reference heating 

rate at stack 43 was Input at nodes 1, 2, and 3. When the analysis was 

performed at stack 37, (Model 38-37-34) the reference heating rate at 
stack 37 was Input at nodes 6 , 7, and 8 , the reference heating rate at 
stack 38 was Input at nodes 11, 12, and 13, and the reference heating 
rate at stack 34 was Input at no es 1, 2, and 3. 

The derived heating ratio q(z)/qief was applied down the walls facing 
the circumferential gap. The same hectinq ratio that was applied under 
node 3 was applied down the wall facing ti'e longitudinal gap for lack of 
better definition. When stack 38 or 39 was being evaluated, the derived 
value of q(z)/qref stacks 38 and 39 was applied down nodes 706, 712, 
718..., and nodes 705, 711, 717..., respectively. At the same time, the 
derived values of q(z)qref at stack 39 were Input to stack 43 (nodes 704, 

710, ?I 6 ...). A different value of heating ratio was not Input to this 

stack because by the time of tests C-188 et al nearly all of the thermo- 
couples near stack 43 had burned through, precluding any possibility of 
obtaining any comparison with temperature In this region. Where stack 37 
was being evaluated (Model 38-37-34) the derived value q(z)/qref was Input 
down stack 38 and the baseline value was input down stack 37 and 34 (In 
addition to being Input down the upper longitudinal gap). 

The reference heating rate at the surface of each of the three thermocouple 
stacks analyzed (37, 38, and 39) for each of the five tests is listed as 
a function of test time in Table II. In order to calculate this heating 
rate the ratio of T/T$*o at each stack was determined by using Figure 3 
which shows T/Ts«o a function of enthalpy for two general locations: 

1) S*-3.0 In. and 2) S»3.0 In., S=7.3 In., and S=7.9 In. The -3.0, 

3.0 and 7.9 In. locations were fr<w) the cold wall model tests of Reference 2, 
and the 7.3 In. location was from the Sllfrax model tests of Reference 1. 
Although the data for three values of S are superimposed on the upper curve 
one line may be faired through all the data points. 

By observation of Figure 4 it may be seen that the value of T/Ts=o is 
approximately the same at all three of the locations in Figure 3. Hence at 
each enthalpy associated with the five cases, the upper curve of Figure 3 
was used to obtain T/Ts»o ^or stack 39. These ratios are tabulated in 
Table I along with the temperatures. Since S =0 at stack 38, T=Ts=o at 
this stack. The value of Ts=o at stack 38 was taken to be Tia (assumed 
constant across the model) and was obtained from the facility printout as 
a function of time. For this temperature (at stack 38) the surface 
heating rate was computed as qref'tJeTw'* (e=.85) In Table II. For stacks 37 
and 39, the ratio T/Ts *0 was multiplied by the temperature at stack 38 to 
obtain T 37 and T 39 (see Table I), and the corresponding value of qref was 
calculated at these temperatures. 
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When stack 37 was analyzed, the thermal model Included stacks 38, 37, and 34. 
To obtain T/T$*0 at stack 34, Figure 4 was used to Interpolate the ratio 
T/T$>0 using a value of S«*5.7 In. Ir: the same manner, when stack 39 was 
analyzed (model containing s acks 38, 39, and 43), Figure 4 was used to 
Interpolate the ratio T/T$»o using a value of S*+5.7 In. No surface 
data was obtained close to these points on any of the three test models 
(Sllfrax, cold wall or RSI) and consequently was not plotted In Figure 4. 

3.0 Results from Analysis 

The results from this Wing Glove model analysis will be discussed In two general 
categories: 1) Best-Fit Heating Ratio and 2) Engineering Verification Heating 

Ratio. Category 1 contained the majority of the emphasis of this study as It 
required considerable Iterations In most cases to obtain the heating ratio that 
would produce the best match of temperature data within schedule constraints. 

3.1 Best Fit Heating Ratio 

The original Orbiter baseline gap heating ratio was adjusted upward at 
each nodal depth down to 3 in. In an attempt to match the temperature 
within 20-30°F at each thermocouple depth for stacks 39 and 38. At 
stack 37, because of the combined low pressures and pressure gradients, 
only the baseline value of J-eating ratio was used. This heating ratio con- 
sistently overpredicted the side wall temperature response at this location, 
possibly because of the presence of nearly separated flow in this region. 

In Figures 5-9 the best-fit heating ratios have been used to correlate 
with temperature data at T/C stack 39 for Tests C-188, 189, 194, 195, and 
190, respectively. Tests C-188, 189, and 195 were for a high stagnation 
pressure (''>22-23 PSF) case. Test C-194 was for a moderate (''>16 PSF) 
pressure case, and Test C-190 was for a low pressure (''>9 PSF) case as seen In 
Table I and Figure 1. The principal time investigated for correlation 
with data was 600 sec. for all cases except C-194 which was 500 sec. In 
addition to these times shown In Figures 5-9, data for additional times 
(100, 200, 300, 400, and 800 sec.) are plotted in some of the figures. 

It may be seen that in most cases the predicted values at 600 sec. 

(500 sec. for C-195) are within 20-30*F of the test data. 

Figures 10-14 present the temperatures at stack 38 using the best-fit 
heating ratios for Tests C-188, 189, 195, 194. and 190, respectively. 

The maximum temperatures In the gap at this stack are some 400-500®F 
lower than at stack 39. This stack had two extra thermocouples to 
correlate with the predictions. For practically all points on all cases 
except the low pressure case (C-190) the predicted tempe''‘ature agreed 
within 20-30“F of the test data for the 600 and 500 sec. times. It was 
determined for C-190 that the best-fit heating ratio curve was actually 
lower than the baseline curve; hence the baseline curve is shown in 
Figure 14 for comparison with the data. 
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The temperatures at stack 37 are shown plotted in Figures 15-19 using the 
zero pressure gradient g(2)/qref Tests C-188, 189, 195, 194, and 190, 
respectively. In this low pressure, moderate pressure gradient region it 
was seen that the baseline heating ratio over predicted the data by as 
much as 200‘’F in some cases. A faired line through the data would be nearly 
parallel to the predicted temperatures for most of the cases for stack 37. 

The uncertainty in the measured surface temperature (from which qref 
was based) in this region could have contributed to some of the lack of 
correlation at this stack. No "best-fit" heating ratio was used for 
this stack because of schedule constraints c^.; the analysis. 

The actual heating ratios used in this best-fit analysis are shown 
in Figure 20 for stacks 38 and 39. There is one curve with the base- 
line value used for stack 38 for Test C-190; this curve was also used 
for all of the stack 37 cases. An additional curve is shown for comparison 
purposes for case C-190 (C-190V) for stack 38 which is lower than the 
baseline value which actually gave a better fit to the data. All of the 
letter designations next to the case number in Figure 20 denote the 
iteration number used to obtain the best fit to the data. 

It may be noted in Figure 20 that a value of PAP/AX is listed after each 
case number. There is an approximate correlation of heating ratio with 
this product, although there may be at least a 10% scatter in determining 
this product. For each stack 39 and 38 there appears to be a definite 
correlation of q(z)/qref with PAP/AX as seen in Figures 21 and 22, 
respectively. For location 38 (lower pressure) the heating ratio drops off 
considerably at low values of PAP/AX while for location 39 (higher pressure) 
the heating ratio is more linear at low values of PAP/AX. 

3.2 Engineering Verification Heating Ratio 

From Figure 20 it may be seen that several of the heating ratio curves 
tend to overlap each other. This is because of the attempt to achieve a 
close fit to experimentally measured temperatures that could be off at 
least 3-5% (as much as 50-80® at Z=0.1 in. and 20-30®F at Z®2.375 in.). 

Also, all of the curves in Figure 20 probably have at least a 10% error in 
the product PAP/AX. The top three curves of stack 39 and 38 (Cases C-188, 
189, and 195) can be represented by a single heating ratio curve that will 
adequately predict temperatures at each thermocouple location. 

In this manner a set of "engineering verification" curves was developed 
at each stack which were approximately parallel to each other and would not 
overlap at each stack. One of these curves was developed for all three 
cases (C-188, 189, and 195). These curves were all designed to at least 
predict or overpredict all thermocouple measurements. Figures 23 and 24 
present these smoothed heating ratio curves at stack 39 and 38, respectively. 
Again the curve for stack 38, Test C-190 is the baseline, or zero pressure 
gradient curve. 
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Figure 25 shows the verification heating ratios plotted as a function of 
PAP/AX for stacks 38 and 39. Only two points are shown for stack 3i' since 
the third point is the baseline value for Case 190 shown in Figure 24. 

This value is actually higher than C-194 (see Figure 24) at lou values 
of Z and then drops to zero at Z»0.40 in. From Figure 25 it may be 
seen that the heating ratios have a higher slope for stack 38 than for 
stack 39. This trend is consistent with the best-fit heating ratios of 
Figures 21 and 22, although there is some scatter in these "non-smoothed" 
heating ratio curves. 

The comparison of predicted temperatures with Wing Glove model data 
using the engineering verification heating ratios are shown in Figures 
26-40. Stacie 39 comparisons are shown in Figures 26-30, stack 38 
comparisons are shown in Figures 31-35, and stack 37 comparisons are 
shown in Figures 36-40 for Cases C-188, 189, 195, 194, and 190, 
respectively. In all of these figures only the comp<irison at 600 sec. 

(500 sec. for Test C-195) is shown. For stack 38 two curves are shown - 
for model 38-38-43 and model 38-37-34 - in order to assess the effects 
of lateral conduction on the temperature predictions. 

A comparison of the temperatures predicted using the engineering 
verification heating ratios with the test data is shown in Figures 
41-46 as a function of test tine of the arc jet. Figures 41, 42, and 
43 present the comparison at stack 39 for Cases C-195, 194, and 190, 
respectively. Figures 44, 45, and 46 show the comparison at stack 38 for 
Cases C-195, 194, and 190, respectively. It may be seen that there is 
better agreement at stack 38 than at stack 39, especially for the 
temperatures at Z=0.1 in. For stack 39 the data was higher tlian the 
predictions during the initial 100 sec. For Test C-190 there was an 
early abort and than a restart wh^ch caused the higher initial temperatures 
(see Figure 54). The in-depth predictions were either close to the data 
or conservative with respect to it for all cases with the exception of 
Test 190 because of the restart. 

4.0 Conclusion 

This document has presented a correlation of predicted temperatures with RSI 
Wing Glove model test data using derived values of in-depth gap heating ratios. 
Best-fit heating ratios were used with a number of iterations to produce a very 
close correlation with the data. Smoothed curves (engineering verification) 
were also developed that were similar in shape to each other and were conservative 
with respect to the data. It was seen that all of the heating ratios correlated 
well with the product of pressure and pressure gradient (low ratios for low 
PAP/AX and high ratios for high PAP/AX for all values of Z). The results of this 
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analysis will be used to make flight predictions of heating in the tile gaps 
at various Orbiter body points and will also be used to compare with recent 
arc-jet Gap/Step (flat plate) model data (References 6 and 7) and with data 
from the NASA/Ames Double Wedge Test Article high AP tests. 
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TABLE II. REFERENCE HEATING RATES USED IN WING GLOVE MODEL THERMAL ANALYSIS 




A. 

Run 

188 




TIME (SEC) 


*-SEC) 


-SEC) 

■'REFjjfWfT'- 

-SEC) 

0 

0.0298 



0.J34 


0.034 


20 

1.78 



3.88 


6.4 


40 

2.40 



5.38 


8.99 


80 

2.27 



5.055 


8.42 


800 

2.27 



5.055 


8.42 







OR; 








or 

I'OOK Q(’ 




B. 

Run 

189 




TIME (SEC) 

qREP3^(BTU/FT^ 

-SEC) 


-SEC) 

qppp (BTU/FT^. 

REF 39 

-SEC) 

0 

0.02792 



0.03193 


0.03193 


30 

2.8046 



6.606 


10.651 


45 

3.0730 



7.288 


11.805 


200 

2.7530 



6.475 


10.446 


800 

2.5539 



5.972 


9.619 



TIME (SEC) 

C. 

qj^PP^^(BTU/FT 2 -SEC) 

Run 190 

qREF (BTU/FT"-SEC) 
38 

qREP (BTU/FT^'-SEC) 
39 

0 

0.0824 

0.135 

0.1576 

20 

2.319 

6.48 

8.71 

60 

3,344 

9.65 

13.09 

80 

3.289 

9.48 

12.85 

H5 

3.457 

10.01 

13.58 

190 

3.457 

10 . Ol 

13.58 

270 

3.752 

10.94 

14.87 

630 

3.813 

11.13 

15.14 

800 

3.937 

11.52 

15.68 
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TABLE II. REFERENCE HEATING RATES USED IN WING GLOVE MODEL THERMAL ANALYSIS 
(CONCLUDED) 


D. Run 194 


TIME (SEC) 


^ref37(btu/ft^-sec) 


(BTU/FT^-SEC) 

38 


q^PP (BTU/FT^-SEC' 
39 


0 

0.029 

0.035 

0.035 

60 

2.35 

6.20 

8.80 

120 

3.25 

8 3 

12.60 

200 

3.49 

9.50 

13.60 

800 

3.49 

9.50 

OF ;‘OOH 



E. Run 195 



TIME (SEC) 


q (btu/ft^-se:) 

Ktl-37 


qpEF (BTU/FT^-SEC) 
38 


q^PP (BTU/FT"-SEC) 
39 


0 

0.029 

0.0344 

0.0376 

10 

1.07 

2.40 

3.58 

20 

2.46 

5.97 

9.22 

40 

3.53 

8.81 

13.76 

80 

3.40 

8.49 

13.25 

530 

3.07 

7.58 

11.78 

1000 

3.07 

7.58 

11.78 
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